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Abstract. 

We discuss the status of thermal model descriptions of particle ratios in 
central nucleus-nucleus collisions at ultra-relativistic energy. An alternative to 
the “Cleymans-Redlich” interpretation of the freeze-out trajectory is given in 
terms of the total baryon density. Emphasis is placed on the relation between the 
chemical equilibration parameters and the QCD phase boundary. Furthermore, 
we trace the essential difference between thermal model analyses of data from 
collisions between elementary particles and from heavy ion collisions as due to a 
transition from local strangeness conservation to percolation of strangeness over 
large volumes, as occurs naturally in a deconfined medium. We also discuss 
predictions of the thermal model for composite particle production. 


The thermal model has been used very successfully to reproduce and (for RHIC 
data even) predict particle ratios measured in central collisions of heavy nuclei at ultra- 
relativistic energy & Si II I- In particular, the predictions for the RHIC data 
are very close to the experimental observations: for a temperature of T=174 MeV 
and a baryon chemical potential of /if, = 46 MeV all available data measured near 
central rapidity are reproduced with a x 2 of 5.7 for 7 effective degrees of freedom 
We note that these calculations are performed under the condition of full strangeness 
equilibration and global strangeness conservation, i.e. the flow of strangeness into and 
out of the central rapidity slice is assumed to cancel. For an update on the inclusion 
of the most recent data see [7j. where it is demonstrated that even the <j>/h~ ratio is 
well described with the same set of parameters. 

These analyses yield a parameter pair (T, fit) for each energy where data have 
been obtained. The current state of affairs is described in Fig. [I]. In this plot 
of temperature vs chemical potential, the location of the phase boundary between 
hadronic matter and the quark-gluon plasma is indicated by the dotted and dashed 
lines, taken from the latest calculations within the framework of lattice QCD || and 
from a bag model equation of state |J. The data points in this figure represent 
chemical freeze-out points determined by a thermal analysis 11 !• & 0 of hadron 
multiplicities measured in central Pb+Pb or Au+Au collisions at various beam 
energies at the SIS, AGS, SPS and RHIC accelerator facilities. The line connecting 
the data points demonstrates, in the framework of the hadron resonance gas model 
used in |j), that chemical freeze-out appears to take place at constant total baryon 
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Figure 1 . The QCD Phase Diagram in the context of the current knowledge 
about relativistic nuclear collisions. For details see text. 


plus anti-baryon density n*, = 0.12/fm 3 . The so determined freeze-out line crosses the 
calculated phase boundary roughly at full SPS energy. 

We note that this interpretation of the freeze-out trajectory differs significantly 
from that of ||, where chemical freeze-out is assumed to take place at constant total 
energy per particle < E > / < N >ps 1 GeV. To us the constant total baryon 
density scenario is more appealing as it implies a physical picture: chemical freeze-out 
takes place at a critical baryon density through baryon-baryon and baryon meson- 
interactions. The absolute value of the critical baryon density n;, = 0.12/fm 3 depends 
of course on the fact that, in J3j, a repulsive hadron-hadron interaction is implemented 
via an excluded volume. 

The fact that the chemical freeze-out points closely approach the calculated QCD 
phase boundary lends strong support to the interpretation @, 0. 0. []. 0 that the 
matter produced in nuclear collisions at SPS and RHIC energies was first thermalized 
in the deconfined quark-gluon plasma phase and subsequently expanded through the 
phase boundary into a thermal gas of (weakly interacting) hadrons. 

Furthermore, as is illustrated in Fig. p], the hadronic matter underneath the 
phase boundary for values of the baryon chemical potential fib < 300 MeV, roughly 
corresponding to full SPS energy and higher, is always dilute: the total baryon density 
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never exceeds that of normal nuclear matter. We, consequently, do not expect a strong 
influence on the chemical freeze-out scenario due to changed hadron masses in the 
medium, if SPS and RHIC energies are considered. Furthermore, even if one includes 
very strong medium effects such as in the chiral model of ]l5|], the chemical freeze-out 
still occurs close to the phase boundary. In the analysis of ||15| the phase boundary is 
at T c = 154 MeV, while the chiral model analysis of the data yields T=144 MeV. 

It is only at lower beam energies, between 5 and 50 GeV/nucleon, that one can 
identify between chemical freeze-out and the phase boundary a phase of dense hadronic 
matter, shown by the hatched area in Fig. |l|. It is near these beam energies, i.e. close 
to 30 GeV/nucleon, where also the relative strangeness content of the hadronic medium 
is maximum jic| . 

This interpretation calls into question the widely accepted belief that the 
enhanced yield of low mass dileptons observed at full SPS energy by the CERES 
collaboration 00 is mainly due to interactions of the p -meson with baryons [jiff . 

At 40 GeV/nucleon beam energy, analysis (2(J of the (still preliminary) hadron 
yields from the NA45, NA49, and NA57 collaborations leads to freeze-out parameters 
of 140 MeV and pb « 400 MeV. This pb value is very close to the value of 
about 420 MeV obtained from the phenomenological parameterization provided by 
equation 1 in [0 at 40 GeV/nucleon. At this energy, the CERES collaboration has 
recently provided first, preliminary results on low-mass di-electrons 00- Again, 
an enhancement relative to the yield expected from standard hadronic sources is 
observed. In fact, the observed enhancement seems to increase from 2.7 ±0.65 at 160 
GeV/nucleon to 4.5 ±1.2 at 40 GeV/nucleon. Whether this is due to the increase, at 
40 compared to 160 GeV/nucleon, of the baryon density between chemical freeze-out 
and the phase boundary is an interesting speculation. However, if Fig. |1] represents 
reality, this increase is only slight, and the real effect should develop at around 20 
GeV/nucleon. Whether the few days of running currently planned at the CERN SPS 
in 2002 at 30 and 20 GeV/nucleon will lead to further insights into this intriguing 
question is unclear at present. 

We note that the planned new GSI accelerator facility [0 is tailored to explore in 
detail the interesting region between 5 and 30 GeV/nucleon, where a dense hadronic 
medium is expected between the QCD phase boundary and chemical freeze-out. 

The above discussed connection between the QCD phase boundary and the 
observed chemical freeze-out points is sometimes called into question |2J, 0 since 
hadron production in e + e _ and pp or pp can also be described in thermal model |f2o| 
yielding a (apparently universal) temperature T e ss 170 MeV. While the fact that T e 
is close to T values determined for heavy ion collisions at top SPS and RHIC energies 
(corresponding to relatively small values of pb ) might indeed reflect the fundamental 
hadronization scale of QCD, we note that there is an essential difference between 
thermal descriptions of central heavy ion collisions and elementary particle reactions. 
As discussed, particle densities and their ratios can be, for heavy ion collisions at 
full AGS energy and higher, well described in the grand-canonical ensemble using the 
particle density n gc . In contrast, for a description of particle production in elementary 
particle collisions, local quantum number conservation needs to be taken into account 
explicitly, leading, for the production of particles with strangeness n to a canonical 
suppression factor 


I n (x i) 
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with which the grand-canonical densities have to be multiplied to obtain the results 
in the canonical ensemble. Here, I n denotes a modified Bessel function of degree 
n and the parameter X\ = 2Vn gc depends explicitly on the volume V over which 
quantum number conservation is enforced |27|. For X\ » 1, —> 1, the 

volume drops out and the grand-canonical limit applicable to ultra-relativistic nuclear 
collisions is recovered. For xi « 1, • This case, appropriate for 

the analysis of e" l “e _ and pp or pp collisions, implies that densities of particles with 
strangeness n contain a quantum number dependent volume V. Hence, the K/-7T ratios 
are proportional to V, the H/7r ratios proportional to V 3 . For the description of 
elementary particle collisions V = V 0 « 7 fm 3 J27| is required, i.e. strangeness is 
conserved within approximately the volume of a nucleon. Moving to an analysis of 

that V = Yq- N part /2 where N part is the number 


27 


heavy ion collisions one finds 
of participants. Already in central Au-Au collisions at SIS energies this “correlation 
volume” V exceeds 1000 fm 3 . The required volume increases further as the beam 
energy is increased, until the grand-canonical limit is reached at full SPS energy. 

This implies that, in central nucleus-nucleus collisions at ultra-relativistic 
energies, strangeness percolates freely over volumes of thousands of fm 3 ! At top SPS 
and RHIC energies it is natural to conclude that the percolation has its origin in the 
quark-gluon phase, lending further strong support to the interpretation above that 
the “coincidence” between experimentally determined chemical freeze-out points and 
the calculated phase boundary implies that a deconfmed phase was produced in such 
collisions. An interesting question is what leads to the percolation at SIS energies, 
where the freeze-out points are far away from where we think the phase boundary 
is. We speculate that, in this case, the percolation takes place in the high density 
hadronic phase indicated by the hatched area in Fig. [|. 

An often overlooked aspect of the thermal model is the possibility to compute also 
the yields of composite particles. For example, the d/p and d/p ratios measured at 
SPS and AGS energies are well reproduced with the same parameters which are used to 
describe baryon and meson ratios @, g, |J. Furthermore, the AGS E864 collaboration 
has recently published Q yields for composite particles (light nuclei) produced in 
central Au-Au collisions at AGS energy near mid-rapidity. In this investigation, an 
exponential decrease of composite particle yield with mass is observed, implying a 
penalty factor of about 48 for each additional nucleon. Extrapolation of the data to 
large transverse momentum values reduces this penalty factor to about 26, principally 
because of transverse flow. In the thermal model, this penalty factor can be easily 
derived. In the relevant Boltzmann approximation, we obtain 


R 


p 


exp 


m± pb 
T 


( 2 ) 


where m is the nucleon mass and the negative sign applies for matter, the positive for 
anti-matter. Small corrections due to the spin degeneracy and the A 3 / 2 term in front 
of the exponential in the Boltzmann formula for particle density are neglected. Using 
the freeze-out parameters T=125 MeV and = 540 MeV appropriate for AGS energy 
J[|, H we calculate R p « 23, in close agreement with the data for the production of light 
nuclei. We also note that the anti-matter yields measured by the E864 collaboration 
p9[ yield penalty factors of about 2-10 5 , again close to the predicted value of 1.3-10 5 . 

This rather quantitative agreement between measured relative yields for 
composite particles and thermal model predictions provides some confidence in the 
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Table 1. Produced number of nonstrange and strange clusters and of strange 
quark matter per central Au+Au collisions at AGS energy, calculated in a thermal 
model for two different temperatures, baryon chemical potential /q, = 0.54 GeV 
and strangeness chemical potential fi s such that overall strangeness is conserved. 
The coalescence model predictions are from Table 2 of |3l| . 
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predictions for yields of exotic objects produced in central nuclear collisions. We 
briefly comment here on the results of a 1995 analysis f30| . 

In this investigation, the production probabilities for exotic strange objects and, 
in particular, for strangelets were computed in the thermal model. The results are 
reproduced in Table [I] for temperatures relevant for beam energies between 10 and 
40 GeV/nucleon. We first note that predictions of the thermal model and, where 
available, the coalescence model of 131 agree (maybe surprisingly) well. Secondly, 
inspection of Table |l] also shows that, in future high statistics experiments which 
will be possible at the planned new GSI facility ^3|, multi-strange objects such 
as I 0 AA^ e should be experimentally accessible with a planned sensitivity of about 
10” 13 per central collision in a years running, should they exist and be produced with 
thermal yields. Investigation of yields of even the lightest conceivable strangelets will 
be difficult, though. 

In summary, we have discussed the surprising success of the thermal model 
to reproduce particle yields in central nucleus-nucleus collisions at ultra-relativistic 
energies. These investigations yield chemical freeze-out parameters which approach 
the phase boundary for energies at and above full SPS energy. The implications of 
this observation are, as discussed above, that the matter produced in nuclear collisions 
at SPS and RHIC energies was most likely first thermalized in a deconfined quark- 
gluon plasma phase and subsequently expanded through the phase boundary into a 
gas of hadrons. In a detailed discussion of various areas of the phase diagram we have 
identified a zone of dense hadronic matter between the expected phase boundary and 
chemical freeze-out which can be studied experimentally with beams in the energy 
region 10 < E; a ;,/A < 40 GeV. Furthermore, we have traced the essential difference 
between thermal model analyses of data from collisions between elementary particles 
and from heavy ion collisions as due to a transition from local strangeness conservation 
to percolation of strangeness over large volumes, as occurs naturally in a deconfined 
medium. Finally, a discussion of composite particle production has emphasized the 
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success of the thermal model also in this area and indicated possibilities to study 
experimentally also rather exotic and interesting multi-strange hadrons. 
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